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ABSTRACT 

The contemporary discoveries of galaxies and gamma ray bursts (GRBs) at high redshift have supplied the 
first direct information on star formation when the universe was only a few hundred million years old. The 
probable origin of long duration GRBs in the deaths of massive stars would link the universal GRB rate to 
the redshift-dependent star formation rate density, although exactly how is currently unknown. As the most 
distant GRBs and star-forming galaxies probe the reionization epoch, the potential rewards of understanding 
the redshift-dependent ratio of the GRB rate to star formation rate are significant and include addressing 
fundamental questions such as incompleteness in rest-frame UV surveys for determining the star formation rate 
at high redshift and time variations in the stellar initial mass function. Using an extensive sample of 1 12 GRBs 
above a fixed luminosity limit drawn from the Second Swift Burst Alert Telescope catalog and accounting 
for uncertainty in their redshift distribution by considering the contribution of "dark" GRBs, we compare the 
cumulative redshift distribution A^(< z) of GRBs with the star formation density /d*(z) measured from UV- 
selected galaxies over < z <4. Strong evolution (e.g., ^'(z) c>c (1 +z)'^ ) is disfavored (Kolmogorov-Smirnov 
test P < 0.07). We show that more modest evolution (e.g., 'I'(z) oc (1 +z)°-^ ) is consistent with the data {P « 0.9) 
and can be readily explained if GRBs occur primarily in low-metallicity galaxies which are proportionally more 
numerous at earlier times. If such trends continue beyond z —4, we find the discovery rate of distant GRBs 
implies a star formation rate density much higher than that inferred from UV-selected galaxies. While some 
previous studies of the GRB-star formation rate connection have concluded that GRB-inferred star formation 
at high redshift would be sufficient to maintain cosmic reionization over 6< z <9 and reproduce the observed 
optical depth of Thomson scattering to the cosmic microwave background, we show that such a star formation 
history would over-predict the observed stellar mass density at z > 4 measured from rest-frame optical surveys. 
The resolution of this important disagreement is currently unclear, and the GRB production rate at early times 
is likely more complex than a simple function of star formation rate and progenitor metallicity. 

Subject headings: gamma rays: bursts - galaxies: evolution 



1. EVITRODUCTION 

The history of star formation in the universe is funda- 
mental for determining the redshift-dependent properties of 
the galaxy population, the production of metals in the in- 
tergalactic medium, and the ionization state of cosmic gas. 
Observational probes of star formation at cosmological dis- 
tances are therefore valuable tools for learning about the bulk 
properties of the univ erse and its contents (for a review, see 
[Robertson et"a^l2010^ ■ 

The current frontier in this field concerns the observa- 
tional determination of the comoving density of star forma- 
tion p*(z) beyond a redshift z —6. This effort is key to un- 
derstanding whether early galaxies were responsible for cos- 
mic reionization as well as when this phase transition in the 
intergalactic medium (IGM) occurred. Good progress has 
been made through efforts to identify rest-frame ultravio- 
let (UV)-selected galaxies in deep optical and near-infrared 
imaging o ut to z ^ 8 (e.g., Jjo uwens et al. 2007, 2010a,b!; 
lOesch et al . 2010; M cLure et al.ll2010l) . Recent spectroscopic 
campaigns have begun to c onfirm th at many of these sources 
do indeed lie at z > 7 (e.g .live et a l. 2006; Stark et al. 20ld; 
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However, the direct detection of z > 6 galaxies spanning the 
full range of luminosities necessary to reliably deduce the in- 
tegrated star formation rate represents a challenge with cur- 
rent facilities. For example, the observationally-inferred /9*(z) 
depends on the adopted value of the poorly-determined lumi- 
nosity f unction faint end slope and the magnitude l imit of the 
survey (iMcLure et al.ll20Tol:lBouwens et alJl2010ah . 

As long-duration Gamma Ray Bursts (GRBs) are thought 
to occur through the demise o f very massive, pos sibly metal- 
poor stars (for a review, see IWooslev & Bloomir2006) . and 
have been observationally connected to broad - line Type Ic 
supernovae (SNe Ic; e.g ., iGalama et alJ Il998t IStanek et"an 
2003; Hior th et all 120031) . the rate of high redshift GRBs 
of this class could provide a valuab le, complernentary 
estimate of Pi,(z) (e.g., Totani 1997; IWiiers et al.l 119981; 
Mao & Mo 1998; Porciani & Madau 2001; Bromm & Loebl 
i2002: ,Chary et aL .2007). In particular, as luminous events, 
they could effectively probe the full extent of the star-forming 
galaxy population including low mass systems undetected in 
the deepest optical and near-infrared surveys. The recent 
discovery of several z > 6 long-duratio n GRBs ( Kawa i et aTl 
2006t iGre iner et al.ll2009t iTanvir et al.ll2009: ^Salvaterraetai] 
2009a; Cu cchiara et al ■1 1201 ih makes this a particularly inter- 
esting avenue to explore. 

The connection between the rate of GRBs (comoving 
Mpc"-* yr"') and the density of star formation (Mq yr"' 
Mpc"-*) can be coarsely parameterized by ^'(z), the ratio of the 
GRB rate to the star formation rate density p^. The ratio \l'(z) 
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could be inferred by comparing the cumulative redshift dis- 
tribution A^(< z) of GRBs with the observed evolution of the 
star-formation rate density over the range of r edshifts 
where is well-measured. Earlier wor k ( Ytik sel et al.l 
Hogs'; 'Kistl eretal.ll2008l I^OOO"; 'W vithe et al.ll2010l) has sug- 
gested that this ratio increases with redshift in the sense that 
GRBs are more frequent for a given star formation rate den- 
sity at earlier times. For instance, if such a redshift depen- 
dence is paramete rized as a simple power law, oc (1 
iKistler et all (|2009() find that ar^ 1.2. Although the origin of 
such a redshift-dependence remains unclear, if extrapolated to 
higher reds hifts the abundance of GRBs at z > 7 (jTanvir et al.l 
12009; .Salvaterra et al.ll2009at ICucchiara et alHIpTlh impUes 
a density of star formation considerably higher (~ 6 - 7 x 
larger) than measured from the abundance of distant UV- 
selected galaxies. In other words, the decline in the GRB rate 
at z > 4 is considerably less rapid than the decrease in the 
star formation rate (e.g., Wanderman & Piran 2010). This dis- 
crepancy has potentially important implications for the reion- 
ization history of the universe and partic ularly the contribu- 
tion from early star -forming galaxies (e.g. . IKistler et al.l2009l; 
lWvitheetal.ll2()T0h . 

Previous determinations of ^(z) at low to intermediate 
redshift (e. g., YuksHetaDEOOS; Kist leret al.1 120081 120091: 
llshida et"ar.li201 1) have not accounted for the uncertainty in 
the GRB redshift distribution A^(< z) that arises from the so- 
called "dark" GRBs defined as those with absent or faint op- 
tical afterglows. Perley et al. (2009) have shown that dark 
GRBs likely span a wide range of redshifts and so they can 
be used to provide a valuable estimate of the uncertainty 
in the overall distribution, corresponding to approximately 
25% in N(< z) at fi xe d z; see their Figure 8 (see also, e.g., 
IGreiner et alllioi ll and lKruhler et al.ll201 lb . In this work we 
use the recently-released Second Swift Burst Alert Telescope 
(BAT) Catalog of GRBs (Sakamoto et al. 2011) to compile a 
comprehensive sample of 1 12 luminous (Liso > lO''' ergs s"') 
long-duration GRBs with known redshifts and upper limits 
dButler et al.1 120071 lIoTot iFvnbo et al.ll2009l: ISakamoto et all 
12011b . accounting for the uncertainty in the overall distri- 
bution by including constra ints provided for a representative 
subsample of dark GRB s (iPerley et all 120091: IGreiner et all 
IMltlKruhler et al.ll20TTl) . 

Our goal is to revisit the redshift-dependent GRB rate to 
star-formation rate ratio 4'(z) with this improved dataset to 
determine carefully whether the ratio evolves and if so on 
what physical basis. A likely driver of such evolution is 
metallicity in the host galaxy population. To explore this 
possibility, we examine whether the evolution of the stel- 
lar mass-metallicity (e.g.. lTremonti et al.l2004tlSavaglio et al.l 
2005; Erb et al. 200^and star formation- stellar mass (e.g., 
Drorv & Alvarezi |2008|) relations, when combined with the 
redshift-dependent stellar mass function (e.g., Drory et al] 
l2005h . can explain the redshift scaling of ^(z). 

Our physical analysis of 'i'iz) over 0< z < 4 enables 
us to calculate the high-redshift (4 < z < 10) star forma- 
tion rate density implied by the presence of distant GRBs 
and to compare this /d*(z) with estimates drawn from rest- 
frame UV galaxy surv eys. Following arguments developed in 
iRobertson et al.l (12010 ) we calculate the impact of this GRB- 
deduced star formation rate density on other indicators of the 
reionization history, including the optical depth of electron 
scattering to the cosmic microwave background and the stel- 
lar mass density that represents the integral of earlier activity. 



We show these additional constraints provide a critically im- 
portant boundary on what is physically plausible in terms of 
the GRB-derived star formation rate at z > 6 in addition to 
the likelihood of the various models we develop to explain 
the redshift dependence of ^(z). We conclude by summariz- 
ing the issues that will need reckoning before the connection 
between the high-redshift GRB and star formation rates is un- 
derstood. 

This paper is organized as follows. We construct a robust 
sample of GRBs above a fixed luminosity limit and discuss 
their redshift distribution and its uncertainties in ^ In O 
we compare the observed cumulative redshift distribution of 
GRBs over < z <4 to the equivalent star formation rate den- 
sity and its evolution. This comparison enables us to consider 
whether and how the GRB rate to star formation rate ratio 

evolves. In we interpret our results in the context of 
a hypothesis where GRBs primarily occur in low-metallicity 
galaxies. We compare predictions of this model with the ob- 
served GRB redshift distribution. In ^ we then use the ob- 
served GRB rate beyond Z2±4 to infer the star formation rate 
density under various assumed forms for ^(z). We discuss the 
ramifications of our results in the context of other constraints 
on cosmic reionization in ^ 

Throughout the paper we have assumed a standard flat 
ACDM cosmology ( fl,,, = 0.3, VIa = 0.7, h = 0.7). 

2. THE REDSHIFT DISTRIBUTION OF LONG DURATION GRBS 

Our goal is to construct a robust sample of GRBs, observed 
down to a fixed luminosity limit, within a redshift range where 
a meaningful comparison can be made with an independently- 
determined star formation rate density, p*(z), from rest-frame 
UV surveys of galaxies. In this way, the connection between 
the GRB rate and the star formation rate density, ^I^, and its 
possible evolution can be studi ed. 

In a manner similar to Kistler et al.l (12009b . we can describe 
the observable number of GRBs within a redshift range Za < 
z < Zfo as 

dV dz 

mza,zb)=K p^izmz)—— (1) 

Jz„ dz 1+z 

where /9+(z) is the global star formation rate density, ^' is the 
number of GRBs per unit star formation rate, dV/dz is the 
redshift-dependent volume element, and the factor 1 /( 1 + z) 
accounts for cosmological time dilation affecting the observed 
GRB rate. The constant K is an factor that accounts for the 
efficiency of the GRB search (e.g., areal coverage, the survey 
flux limit, etc.), but its value is unimportant for our purposes^ 
We can remove the dependence on K by simply constructing 
the cumulative redshift of GRBs over the redshift range < 
z < Zmax, normalized to A^(0,Zmax), as 

jV«zkmax)= (2) 
A^(0,Zmax) 

The product /9*(z)^'(z) then sensibly determines the cumula- 
tive redshift distribution of GRBs. Given this normalization, 
we will also sometimes refer to the redshift-dependent frac- 
tion of star formation that can produce GRBs iJj(z), which is 

* A clear concern is the possible multiplicative degeneracy between K and 
^(7) in Equation[T] and that possible redshift dependence in, e.g., the follow- 
up efficiency could mimic an evolution of ^'(z). We note that Dai 1 2009) have 
found that the cumulative peak photon flux distribution of Swift GRBs with 
and without spectroscopic redshifts are similar. This result suggests that the 
follow-up efficiency does not strongly depend on redshift. 



related to the number of GRBs per unit star formation rate as 
^'(z) = ^'oV'(z)- The constant ^I^o then encodes the number of 
GRBs formed per unit mass of stars. The value of ^I^o cannot 
be determined independently of the unknown K in Equation 
[T] but the time or redshift when 'So is defined does matter 
for a model where ifjiz) is calculated directly (see ^ below) 
rather than averaged over some redshift interval (c.f.. Equa- 
tion |2l). In such a case, we take to be defined relative to 
when ijj{z) = 1. 

2.1. GRB Catalog 

To evaluate the possibility of redshift-dependence in the 
GRB to star formation rate density ^'(z) through Equation 
[U we require an observational sample to construct the cu- 
mulative redshift distribution calculated by Equation |2l The 
primary requirements are a well-understood completeness 
in the redshift determinations above some GRB isotropic- 
equivalent luminosity. The distribution A^(< zjzmax) can be 
determined from catalogs of GRBs monitored with gamma- 
ray satellites (e.g.. Swift; Gehrels et al. 2004) and followed 
up from the ground (e.g., Fynbo_et al. 2009). For our base 
catalog of GRBs, we take the union of the s amples pre- 
sented inlButler et al J(l2007h.lPerleY et al.' ('2009').'Butl er et alj 
(1201 Oh. IsTk amoto etal.1 (1201 Ih TIOreiner et al. (201 ij), and 
iRriihler etalj(i201 li) . We include only GRBs occurring before 
the end of the Second Swift BAT GRB Catalog and prefer the 
most recent redshifts for GRBs where the samples disagree. 
This union provides a sample of 164 GRBs with known red- 
shifts and redshift upper limits, but two GRBs (071 1 12C and 
060505) have incomplete fluence or burst duration measures 
and are discarded. The remaining 162 long duration GRBs 
with redshifts or redshift limits serve as our base GRB cata- 
log. 

To account for the incom pleteness owing to the flux limit of 
the Swift survey, we follow lKistler et all (12009 ) and construct 
a subsample of with isotropic-equivalent luminosities Liso > 
10^' ergs s"'. The luminosity is computed as 



(3) 



f90/(l+z) 



where Eiso is the isotropic-equivalent energy, fgo is the burst 
duration containing from 5% to 95% of the total fluence, and 
the f actor of ( 1 + f) acco unts for cosmological time dilation 
(e.g.. iKistler et al.lf2009h . For isot ropic-equivalent energies 
129 GRB £iso values are taken from lButler elaH (120071 l20Toh 
and 4 recen t values (0 60908, 090926B, 091018, and 091029) 
from Sakamoto et al.l (1201 lb . Additionally, the isotropic en- 
ergy values for 21 further GRB s are computed from the flu- 
ences reported by iButler et"an (12007. 2010). The isotropic 
energies from the 8 remaining GRBs (060512, 0908 14A, 
090904B, 090927, 091020, 091024, 091127, and 091208B) 
are calculated from the 15-150 keV fluences reported by 
[Sakamoto et al. (201 1), but are possibly lower limits given the 
typical energy range of 1 - lOOOOkeV for de fining isotropic- 
equivalent quantities (e.g.. lAmati et al.|l2002 [). The burst du - 
rations fgo are taken from Table 1 of Sakam oFo et al.l (1201 II) . 
except for GRBs 05 0820A, 060218, and 090529A taken from 
( IButler et all [20071 [2010). For our catalog of 162 GRBs 
with redshifts and redshift limits, this culling provides 112 
GRBs with isotropic-equivalent luminosities Liso > 10''' ergs 
s"' for our analysis. The redshifts and limits, isotropic 
equivalent energies and luminosities, and burst durations 
of the full sample of 162 GRBs compiled from the union 



of the IButler et al.l (120071) iPerlev et al.l (!2009h, IButler et al.l 
(2010), Sakamoto S^f UOl ll) . iGreiner et all (1201 ll) and 
IRriihler et al.l (|201 li) catalogs are provided for convenience in 
Table [TJin the Appendix. 

Since we will use the cumulative redshift distribution A^(< 
z) of this sample as the basis for our analysis, it is important 
to consider its uncertainties. While the Swift catalogs provide 
a valuable compilation of gamma-ray detections, the redshift 
determinations are clearly influenced by their optical observ- 
ability. The phenomenon of so-called "dark" GRBs with sup- 
pressed optical counterparts could influence whether the ob- 
served j y(< z) is rep r esent ative of that for all long-duration 
GRBs. iPerlev et"an (12009 ) have considered this important 
issue by attempting to constrain the redshift distribution of 
dark GRBs through deep searches that successfully located 
faint o ptical and near-infrared counterparts. The Perle y et aP 
(l2009h work provides us with 2 redshifts and 9 redshift up- 
per Umits fo r a sub sam ple of dark GRBs i n our catalog. 
IGreiner et al] (120111) and iKrlihler et all (120111) have pursued 
this effort in parallel, and have provided 3 additional red- 
shifts and 1 redshift upper limit for dark GRBs in our catalog. 
We assume the subsamples of dark GRBs with redshift upper 
limi ts presented b y Perley et al. (2009), Greiner et al. (201 1), 
and [kriihler et al.l (|201 ll) are representative of that class, and 
therefore optionally incorporate those limits to characterize 
the effects of possible incompleteness of the Swift sample with 
firm redshift determinations. We also note that while th e lumi- 
nosity limit for our sample was chosen to match Kistle retaD 
(2009), at redshifts above z > 4 the Swift sample is incom- 
plete for this limit. However, fully accounting for this incom- 
pleteness would only increase the relative number of GRBs at 
high-redshifts. As our following results show, our sample's 
luminosity limit is therefore conservative. 

3. COMPARING GRB RATES TO THE COSMIC STAR FORMATION 
HISTORY 

As we have yet to develop physical intuition into the con- 
nection between the rate of GRBs and /9*(z), we will begin 
by an empirical comparison of the cumulative GRB red shift 
distribution A^(< z|zmax) constructed as described in ^2.1l with 
the cumulative redshift distribution that would be expected 
given the observed star formation rate density /i*(z) from rest- 
frame UV-selected sample and various forms for the redshift- 
dependent ratio 5'(z). 

For p*(z) we use the results from iHopkins & Beacom[ 
(l2006h who gathered and standardized measures of the 
star formation rate density from the Hopkins (200^) 
compilation a nd observatio ns by Wolf et al. (2003|, 
'Bouwens_et_al.| (2003a b) . Bunker et al.l (2004), Ouchi et aU 
(2004), Arno uts et alij (120051), [l TfIoc' h et al. (2005E 
Perez-Gonzalez et al. " ('2005'), 'Schiminovich et al. (2 005|) . 
Bouwens & Illingworth (2006), Hanish et al. (2006), and 
Thom pson et al.i (i2006i) . iHopki ns & Beacom (2006) provide 
a piecewise-linear "Modified Salpeter A IMF" model in their 
Table 2 that provides a good statistical fit to the available 
star formation density data. We limit our use of their fit 
to z < 4 where the data is optimal. We note that p*(z) is 
relatively flat [p*(z) oc (1 +z)""-^'^] in the interval 1 < z < 4 
where most of the GRBs with Liso > 10^' ergs s"' occur. 
This scaling means the use of Equation [2| is particularly 
accurate as the factor of ^ 2 in the normal ization of p<,(z) 
allowe d by the 3 - tr uncertainty in the Ho pkins & Beacoml 
( 120061) fit is circumvented (see also ^ below). We note 
that using the [Hopkins & Beacom (.2006.) results for the 
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iBaldrv & Glazebrookl (l2003h IMF instead of the Salpeter A 
IMF model (see, e.g. iHopkins & Beacoinl 120081) has little 
effect on our conclusions. 

Figure [T] shows the cumulative redshift distribution of ob- 
served GRBs (black histogram), normalized over the redshift 
range < z < 4. The gray-shaded region shows how the distri- 
bution shifts in the limiting cases of all dark GR Bs occurring 
at z = or the upper redshift limits determined by Perlev et aLl 
(l2009l) . iGreiner et al.1 (l20Tll) . and iKruhleretalJ OOllI) . Fig- 
ure[T]compares the observed GRB cumulative redshift distri- 
bution A^(< zjzmax = 4) for the case of three models for the 
redshift evolution of the GRB rate to SFR ratio 5'(z). 

If the quantity 'i'(z) ^ is constant (red line), the cumulative 
redshift distribution of GRBs increases rapidly at z 2-3 in 
sync with the star formation rate density. If instead the GRB 
rate to SFR ratio evolves strongly with redshift over the epoch 
z < 4 as, for instance, ^'(z) oc (1 + z)'^ (orange line), then the 
GRB rate is shifted to higher redshifts and the cumulative dis- 
tribution increases rapidly at z > 3. A weaker redshift evolu- 
tion [^(z) oc (1 + z)°'^] better reproduces the cumulative GRB 
rate density (blue line). Given the relatively small sample, 
the data appears roughly consistent with each of these ^'(z) 
redshift-scalings. 

3.1. Statistical Tests 

Given the integral distribution of observed GRB redshifts, 
and a parameterized model for predicting this distribution, 
we can perform statistical tests to check for consistency be- 
tween the observed and model distributions, and calculate 
confidence intervals for the parameters of the model given the 
data. The null hypothesis that the observed GRB redshifts are 
consistent with a model distribution can be evaluated by the 
one-sample Kolmogorov-Smirnov (KS) test, which calculates 
a "f -value" that corresponds to one minus the probability that 
the null hypothesis can be rejected. We employ this statistical 
test to determine a plausible range of values for the paramter 
a in a model where 5'(z) oc (1 +z)". We only use the KS test 
to conservatively evaluate the relative agreement between the 
observed and model GRB redshift distributions. More for- 
mally, we can calculate the likelihood function for the param- 
eter a given the observed data and assumed parameterized 
model using a Bayesian technique. For the case of inde- 
pendent samples z = [z,] from a redshift probability distribu- 
tion p(z\a) (e.g., the integrand of Equation[TJ, the likelihood 
function of a can be calculated as 



N-\ 



£.ia) = Y[p(zi\a). 



(4) 



(=0 



The posterior probability density of a given the observed data 
Zi and parameterized model ^ can then be computed as 



p(a|z,*): 



p(a\Wa) 



(5) 



where p(a 1 5*) is the prior probability of the parameter a given 
the model ^P. We take this prior to be flat over a wide range of 
a such that it does not affect the shape of pia\z, ^f). The con- 
fidence interval corresponding to a given Gaussian equivalent 
significance can then be calculated by integrating Equation 
|5] about the peak likelihood, and whenever an effective a is 
quoted it refers to the Bayesian confidence region. 

Figure |2] shows the KS-test probability assuming 'I'(z) oc 
(1 +z)" for -1 < a < 2.5. Comparing only with GRBs with 
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Figure 1. Cumulative redshift distribution of long-duration gamma-ray bursts (GRBs) 
at z < 4 determined from the second Swift BAT catalog with isotropic equivalent lumi- 
nosities Liso > 10^' ergs s"' (black histogram, see text for source of data). Incorporating 
the redshift constraints of dark GRBs shifts the distribution over the range indicated by 
the gray shaded area, and provides a estimate of the intrinsic uncertainty in the obser- 
vational determination of the distribution. Three simple power-law parameterizations of 
the ratio '^{z) between the GRB and star formation rate densities are shown: a constant 
ratio with redshift (red line), ^(z) oc (1 +z)°-' (blue line), and *(z) oc (1 +zf '^. 

spectroscopic redshifts (the solid black histogram in Figure 
[T]), we find that the region where KS P > 0.05 contains power- 
law indices of -0.2 a < 1-5. The peak probability occurs 
for a « 0.5. Computing the posterior probability /?(a|z,\I'), 
we find a constant 4* (i.e., no evolution) is marginally al- 
lowed at the 2- a level. Including the dark GRB redshift 
constraints shifts the probability curve. If all dark GRBs are 
local, then the KS test region P > 0.05 contains power-law 
indices -0.7 < a < 1 .0 (with a median probability near a = 
corresponding to no evolution). Instead, if all dark GRBs are 
at their maximum possible redshift then within 2 -a power- 
law indices 0.3 < a < 1.7 are favored. While the constraining 
power of the current sample is not particularly stringent, fur- 
ther monitoring of GRBs should soon definitively rule-out a 
'I'(z) that is constant or declines with redshift. We discuss the 
implications of this constraint further in ^ 

In summary, we find that the distribution of GRBs with 
spectroscopic redshifts is consistent with only moderate vari- 
ation of ^'(z) over < z ^ 4 and that there is overall only 
weak evidence for evo lution (^ 2g confid ence). Compared to 
previous studies (e.g.. lKistler et alj |2009) the results are con- 
sistent at the ^ 2cr-level, but we infer a weaker redshift depen- 
dence owing to the fractionally increased number of GRBs at 
z < 2 in our compiled GRB sample. The additional uncer- 
tainty arising from including constraints from dark GRBs is 
important to include. If dark GRBs occur at their maximum 
allowed redshifts the distribution is more heavily weighted to- 
wards higher redshifts, and more strongly indicates a possible 
redshift dependence to ^'(z). 

Clearly, to make progress it would be helpful to have a 
physical basis for evolution in 'I'(z). The most Ukely hy- 
pothesis links the GRB production rate with the metallicity 
of the underlying stellar population (e.g.. lWooslev & BloomI 
120061) . Since the typical metallicity of galaxies at fixed stel- 
lar mass dechnes with redshift (e. g., iTremonti et al.l I2004t 
iSavaglio et al.ll2005l : lErb et al.ll2006l) this hypothesis provides 
a natural basis for an increase in \l/(z). 
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+ Dark GRBs @ z=0 
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Figure 2. Results of a Kolmogorov-Smimov test comparing the GRB cumulative 
redshift distribution N(< z\Zmii\ = 4) with model distributions calculated assuming 
the Hopkins & Beacom 1 2006) model for the star formation rate density p<,(z) with a 
redshift-dependent GRB rate to SFR ratio '^(z) oc (1 +^)". Shown is the probabil- 
ity P of consistency as a function of the power-law index a. Comparing with sample 
with spectroscopic redshifts (red line), the P > 0.05 interval is —0.2 < a < 1.5 with 
a peak probability near a = 0.5. The Gaussian equivalent 2-cr confidence interval is 
0.025 < Q < 1.5 (see text). If the dark GRB sample with redshift upper limits is as- 
sumed to be local {z ^ 0), the P > 0.05 interval shifts to -0.7 < a < 1 with a median 
near a ^ 0. If the dark GRBs lie at their maximum possible redshift the interval moves 
to < a < 1.8 with a peak near a ^ 0.8. 



4. MODELING THE REDSHIFT-DEPENDENT GRB TO SFR RATIO 
WITH METALLICITY EVOLUTION 

We now consider a physical scenario where the cosmic 
GRB rate is enhanced in low metallicity galaxies, and de- 
velop a method for calculating the fraction of the star forma- 
tion rate density occurring below a characteristic metallicity 
12-l-log[0/H]crit above which GRB production is suppressed. 
Using the same statistical comparisons employed in ^ we 
can evaluate whether the observed GRB rate is consistent with 
such a picture and, if so, what values of 12-l-log[0/H]crit are 
favored. 

A variety of theor etical pictures for the origin of long 
duration GRBs (e.g.. IWoos"ievlll993h suggest that the GRB 
rate may be connected to the metallicity of their progen- 
itor stellar population ( for reviews, see iMeszarosI 120061: 
iWooslev & Bloo in 2006). GRBs may require the retention 
of significant angular momentum after collapse, which limits 
the amount of mass loss prior to collapse (e.g., Frv er et alj 
119991: iMacFadven et all 1200 Ih . The lower opacity and 
mass loss rates of low-metallicity stars eases this require- 
ment and provides a credible rationale for preferring low- 
metal l icity GRB hosts fe irschi et al. 2005; Yoon & Lang^ 
I200I IWooslev & Hegea T2l 



12006 ). Moreover, host galaxies 
of long duration GRBs are often observed to be metal 
poor and ha ve low luminos i ties, both in the nearb y uni- 



verse 



iModiaz et al. 



^chaska et af '2004'; 'Sollerman et al.' '2005| 
2006; Stanek et al. 2006; Wiersema et al. 200 
and at cosmological distances (e.g., Fruchter et al. 1999, 
2006; Le Floc'h et al. 2003; Fynbo et al. 2003; Savaglio et al. 
20091) . 



The quantitative details of the metall icity connectio n are 
still debated (for a good discussion, see lModiazir201 lb . For 
instance, by exploiting the connection between SNe Ic dis- 
playing broad fines and GRBs (e.g., .Galama et al.. .1998; 



iStaneket al.ll200l iHiortfiet al.ll2003h . IModiaz et all (l2008h 
fou nd tfiat only low-metalficity (12 -I- log[0/H] < 8.7 on 
tfie 'Kobulnicky & Kewley' '2004' scale, converted following 
iKewley & Ellison 2008) galaxies have GRBs that fade into 
SNe Ic, whereas host galaxies with metallicities above this 
ceiling have GRB-free SNe Ic. However, there are regions 
within GRB ho st galaxies known with higher metallicities 
dLevesque et al.''2010b), and the radio-relativistic SN 2009bb 
([Soderberg et al. 2010 1) occurred in a high - metalficity region 
of its host galaxy (Levesaue et al.l l2010d) . iHan et alJ (120101) 
and Levesaue et al. (2010a) find that GRB host galaxies lie 
below the ITremo nti et al. (2004) stellar mass-metallicity re- 
lation, whereas iSavaglio et al. (2009.) report that the prop- 
erties of GRB ho sts do not clear l y diff er from normal star- 
forming galaxies. iMannucci et al.l (1201 ll) find that GRB hosts 
are offset from the mass-metallicity relation owing to their 
their higher star formation rates than average at fixed stel- 
lar mass, but remain within the "fundamental metallicity rela- 
tion" between star form ation rate, stellar mass, and metallicity 
(IMannucci et al.ll2010i) . 

Given this controversy, we considered it appropriate to con- 
struct a model to calculate 5'(z) in the case where GRBs pref- 
erentially occur in host galaxies with metallicitie s below some 
characteri stic l2-l-log[0/H1cri t- iKocevski et al.l (|2009) com- 
bined the IModiaz et al.l (120081) empirical host galaxy metal- 
licity ceiling for GRBs with t he redshift evolution of the stel- 
lar mass-metallicity relation ('Savaglio et al."2005'), the stellar 
mass- star formation rate relation evolution (Drorv & Alvarez 



2008), and the stellar mass function evolution (DroryetaL 
2005) to calculate the redshift-dependent cha racteristic mass 
of GRB host galaxies. Below, we extend the iKocevski et al.) 
(l2009i) formalism to model the redshift dependence of \['(z) 
and allow for a variable metallicity ceiling 12-l-log[0/H]crit- 
We note that a similar approach has been adopted by 
Langer &Norman (2006), Salvaterra & Chincarini (2007^ 
Salvat erraetalJ (2009b), Butler et al. (2010) and Virgil iet all 
(201 L) . Analyses that used the Langer & Norman (20061) cal- 
culation utilized the fraction of stellar mass at metallicities be- 
low 12-l-log[0/H]crit to determine the redshift-dependence of 
5'(z), whereas we prefer to extend the Kocevski et al. (2009]) 
model to find '^{z) from the fraction of star formation occur- 
ring at metallicities below 12-l-log[0/H]crit. 

4.1. A Model for '^{z)from Metallicity Evolution 

Given the relation between stellar mass and metallic- 
ity (e.g., ITremonti et al.l 120041) a metallicity ceiling 12-1- 
log[0/H]crit for GRB host galaxies would imply a critical 
galaxy stellar mass M^^nt above which GRB production is 
suppressed. The redshift-dependent stellar-mass metallicity 
relation can be parameterize d as ( Savaglio et aL.2005„ on the 
iKobulnicky & Kewleyll2004i scale) 

1 2 -I- log[0 /H] = -7 .5903 + 2.53 1 5 logM,, 
-0.09649 log^M,, 
+5.17331ogf„-0.39441og2f„ 
-0.403 logf„logM,, (6) 

where f„ is the age of the universe at redshift z in Gyr and M,t 
is the galaxy stellar mass in solar masses. Equation|6]then can 
be used to connect a given metallicity ceiling 12-l-log[0/H]crit 
to a critical galaxy mass M^,,ail{z) for GRB production. 

We now introduce the fraction of star formation occurring 
in galaxies with metallicities lower than 12 -I- log[0/H]crit, 
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which can be expressed as 



Jo 



SFR(M,z)$*(M,z)dM 



/o°°SFR(M,z)$*(M,z)dM 



(7) 



where SFR(M,z) is the star formation rate - stellar mass 
relation and $j,(M,z) i s the galaxy stellar mass function. 
iDrory & AlvarezI (l2008h parameterize the observed star for- 
mation rate-stellar mass relation as 



SFR(M,z) = SFRo 



M 
Mo 



exp 



M 



(8) 



where /3 k, 0.5 and the parameters SFRq and Mq evolve with 
redshift as 



SFRo(z) = 3.01(1 +zy-'"MQ yr' 
Mo(z) = 2.7x 10'°(l+zf Mq. 



(9) 
(10) 



The stellar mas s function also evolv e s with redshift (e.g., 
IDrory et al.ll2005l) . Taking a lSchechted (Il976h model for the 
galaxy stellar mass 



$*(M,z) = 



M 



Ml) ^""^y Ml 



M 



dM 

Ml 



T^, (11) 



iDrorv & Alvarez" (2008) model the redshift-dependence of the 
observed galaxy stellar mass function through the parameters 



(t>Jx) ~ 0.003(1 Mpc-^dex"' 
log[Mi/M0](z)« 11. 35-0.22 ln(l+z) 
7(z)«-1.3. 



(12) 
(13) 
(14) 

By combining Equations l6]and[8lfT4l the redshift-dependent 
ratio of the GRB rate to the star formation rate ^'(z) can be 
estimated by evaluating Equation 5. With ^'(z) in hand, a 
model for the cumulative redshift distribution of GRBs can 
be calculated using Equations [T] and |2] once a model of the 
star formation rate density p*(z) is adopted. For a critical host 
galaxy metallicity 12 + log[0/H]ci-it above which GRBs are 
suppressed, we can use the GRB data to inform us as to what 
metallicity ceilings are plausible. The redshift-dependence 
of '^(z) will clearly depend on the value of 12 + log[0/H]ciit 
since the fraction of star formation occurring at metallici- 
ties below 12 + log[0/H]ciit will vary with redshift owing to 
the evolution of the mass-metallicity relation, the mass-star 
formation rate relation, and the stellar mass function. For 
simplicity in this model, GRBs are prevented from occuring 
above 12+log[0/H]crit, but as noted above in ^GRBs do oc- 
cur in metal rich host galaxies. The suppression in metal rich 
galaxies should therefore be understood to be an incomplete, 
coarse model to indicate the potential preference for GRBs to 
occur in metal-poor systems. 

Figure [3] shows the cumulative GRB redshift distribution 
A^(< zkmax = 4) resulting from Eq uation [7] with the adoption 
of the Ho pkins & BeacomI (120061) star formation rate den- 
sity model, for three choices for the critical metallicity 12 + 
log[0/H]ci-it. One model adopts a large 12 + log[0/H](;rit = 
9 (red U ne), sim ilar to the metallicity of the host of GRB 
020819 (iLevesq ue et al. 2010b). In this case, essentially 
all star formation occurs in hosts with metallicities below 
12 + log[0/H](;rit and is roughly constant with redshift as 
the GRB rate and star formation rate density trace one another. 
The large 12 + log[0/H]crit model therefore closely tracks the 
a = model shown in Figure[T] An intermediate model adopts 



5^ 0.6 




GRBs with known z 

including dark GRBs 

if GRBs occur at: 
12+log O/H < 9 
12+log O/H < 8.7 



0.0 1 



Figure 3. Cumulative redshift distribution of gamma-ray bursts at < 4 and 
metallicity-evolution models for the GRB rate to star formation rate ratio '^(z). The 
black histogram and gray area indicate the cumulative distribution of GRBs with firm 
redshifts and the uncertainty owing to dark GRBs {see Figure [Tlcaption). Also shown 
are three models of the ratio '^(z) where GRB production is suppressed in galaxies with 
metallicities above a ceiling 12 + log[0/H]ciii. determined by model presented in ^ 
a high 12 + log[0/H]„ii = 9 (red line), the 12 + log[0/H]crii Ri 8.7 ceiling found by 
Modjaz et al. (2008, blue line), and a illustrative low ceihng 12 + log[0/H]cri, = 8 (or- 
ange line). The inset shows the redshift-dependence ip{z) corresponding to the model 
with 12-t-log[0/H]crit ^ 8.7, along with a parameterized fit (see Equation ! 151 . 



the value of 12-l-log[0/H]crit « 8.7 from lModiaz et all (120081) . 
shown as the blue line in the Figure [3] Star format ion oc- 
curring in galaxies with metallicities below the Modia zet alj 
(2008) 12-l-log[0/H]crit tracks the GRB rate with surprising 
fidelity, and for convenience we show the corresponding '^{z) 
in the Figure [3] inset and provide a parameterized fit to this 
^'(z) model as 



*fi,(z) = 0.5454-1- (1 -0.5454) x [erf(0.324675z)] 



1.45 



(15) 



(dashed black line in Figure [3] inset) that recovers the com- 
puted 7/)(z) to within 1% atO<z< 10. Normalized over the 
redshift range < z < 4, this intermediate 12-l-log[0/H]crit 
model produces a cumulative redshift distribution similar to 
the ^'(z) oc (1 +z)"'^ model discussed in ^ Third, we show 
the effects of a low value of 12 -I- log[0/H]crit = 8. While 
the fraction of star formation occurring in systems with 12-1- 
log[0/H] < 8 is much smaller than 12-l-log[0/H] « 8.7, nor- 
malized over the redshift range < z < 4 the redshift depen- 
dence of the two models are similar Sensibly, the low 12-1- 
log[0/H]ciit model evolves with a somewhat stronger redshift 
dependence as the epoch at where the characteristic stellar 
mass in the stellar mass function reaches 12-l-log[0/H]ci-it is 
pushed to higher redshift. This low 12-l-log[0/H]crit model 
only serves as a strawman to illustrate the calculated redshift 
dependence; GRBs are observed to occur at larger metalli ci- 
ties (see, e.g.. ILevesque et alteOlOatlMannucci et al.ll201 lb . 

As in ^ we can formalize this comparison using a one- 
sample KS test. Figure |4] shows the KS test probability as 
a function of 12 -I- log[0/H1t^rit in ter ms of 12 -I- log[0/H] 
on the Kobulnicky & Kewlev (12004 scale for 8 < 12-1- 
log[0/H]crit < 9.2". We find that a cumulative GRB red- 
shift distribution produced by a 12-l-log[0/H](;rit ^ 8.85 pro- 
duces a redshift-evolution in ^'(z) that is adequately consis- 
tent (P > 0.05) with the firm GRB sample (red line), or the 
sample enlarged by dark GRBs at their redshift limits (or- 



7 



GRBN(z): 

Known Redshifts Only 
+ Dark GRBs @ z=0 
GRBs @ z= 




8.4 8.6 8.8 

12 + logO/H [KK04] 

Figure 4. Kolmogorov-Smimov test results which compare the observed GRB cu- 
mulative redshift distribution with the calculated distribution A^(< z|2max = 4) assuming 
the Hopkins & Beacom i2006) model for the star formation rate density p^,(z) with a 
redshift-dependent GRB rate to SFR ratio '^(z) given by Equation^ Shown is the prob- 
ability P of consistency as a function of the critical metallicity 12-t-log[0/H]i;Hi (on the 
scale of Kobulnicky & Kewley 2004) above which GRBs are suppressed. For the sam- 
ple with firm redshifts (red line) all critical metallicities produce KS P > 0.05, meaning 
the data are consistent with no metallicity ceiling (note that at 12 -t- log[0/H]crii ^ 9 
the models are equivalent to a = in Figure[2). However, the peak probability is very 
similar to the metallicity ceiling claimed by Modjaz et al. 1 2008) (dotted line) and the 
consistency of the model and data rapidly increases for 12-t-log[0/H]LTii < 8.9. Includ- 
ing the dark GRB sample assuming they lie at their maximum possible redshift shifts 
the P > 0.05 interval to 12-t-log[0/H]i;Hi < 8.9 (formally, this interval equates to a 2-cr 
Gaussian equivalent confidence; see text). If dark GRBs were all local (:: ^ 0, blue line), 
no clear value for 12 + log[0/H][:rii is favored. 

ange line). Using Equations [T] and [T] to construct the like- 
lihood function >C(12-l-log[0/H]crit) and the posterior distri- 
bution /9(12-l-log[0/H]crit|z,*), we find that the firm GRB 
redshift sample is formally 2- a consistent with no metallic- 
ity ceiling. Similarly, we find that including the dark GRB 
sample at their redshift limits has a 2- a confidence interval 
of 12-l-log[0/H]ciit < 8.9. We note that in both these models 
the maximum probability occurs for a metallicity si milar to 
the cei ling 12-l-log[0/H]crit ~ 8.7 suggested by Modia Fet alj 
(|2008, Figure |4] dotted line). While the cumulative GRB 
distribution prefers a 12-l-log[0/H]crit S 8-7, lower critical 
metallicities display similar consistency. We note that these 
low critical metallicities can in principle be differentiated 
based on absolute comparisons of the GRB rate, rather than 
through normalized cumulative distribution functions. Fi- 
nally, if dark GRBs are local (z ^ 0) phenomena (blue line), 
then the observed GRB distribution is not very constraining 
with a large range in 12-l-log[0/H]ciit displaying similar con- 
sistency with the data. 

In summary, it is relatively straightforward to construct 
physically-plausible models for ^'(z) within the sample un- 
certainties, based on a metallicity ceiling for GRB produc- 
tion. As the sample sizes grow, there is every prospect of 
securing valuable constraints on such models. As a note, we 
caution that the recent observations by lSavaglio et alJ (1201 Ih 
of the afterglow of GRB 090323 at redshift z ~ 4 shows ev- 
idence for two DL As with supersolar metallicities, and that 
ICenko et al.l (|201 ll) suggest at least one of these DLAs may 
be associated with the GRB host galaxy. If these observations 
are confirmed, they could pose difficulty for models of the 
GRB rate with a low 12-l-log[0/H]ciit. 



5. THE HIGH REDSHIFT STAR FORMATION RATE DENSITY 
DERIVED FROM DISTANT GRBS 

An exciting development in the past few years has been 
the discovery and verification of the first sample of long- 
duration GRBs beyond redshifts z c±6 (e.g ., Kawai et al. 200m 
Tanvir etal.1 120091: ISalvaterra et"aLl l2009a; .Cucchiai-a et all 
201 Ih . We now utilize our understanding of the ratio of the 
GRB rate to the star formation rate density and its pos- 
sible redshift dependence to interpret these data. Of partic- 
ular interest is how such GRB-derived estimates of the star 
formation rate density jO*(z) compare with those determined 
from UV-selected galaxy samples, as these quantities hold im- 
plications for the timing of cosmic reionization and whether 
the density of star-forming galaxies alone provides sufficient 
energetic radiation to reionize the IGM. In the following we 
will f ollow clos ely the procedures and arguments developed 
in Rob ertson et a l. (2010). 

In the manner of .Yiiksel et al.l (120081) . we can estimate the 
star formation rate density as 



(P*)(Zl,Z2)^ 



Ngrb(zuZ2) I' Mzmz)%^^ 

a^grb(i,4) J;:mz)f^ 



(16) 



where AfeRB(zi,Z2) is the observed number of GRBs at red- 
shifts zi < z < Z2- Panel a of Figure |5] shows the resulting 
comparison. The star fo rmation rate density d etermined from 
UV-selected galaxies (McLure et al. 2010; Bouwens et alJ 
|2007, 2010b; Oesch et al. 2010; Schiminovich et al. 20051 
iReddv & Ste idel 200 9D increases to z 3 and then declines 
to high redshift CBouwens et al.ll2010al gray points with er- 
ror bars). Two param eterized star formation histories from 
[Robertson et al.l (12010 1) consistent with the data are shown for 
illustration (green hatched area). Star formation rate densities 
estimated from the high-redshift (z > 4) GRB rate, as calcu- 
lated by Equation [T6l and three models from Figures [T] and [3] 
are shown as red (a = 0, equivalent to 12-l-log[0/H]crit ^ 9), 
blue (12-l-log[0/H]crit ~ 8.7, equivalent to a w 0.5 at z < 4 
and a w at z > 4), and orange points (^'(z) oc (1 +z)''^, 
stronge r evo lution than any metallicity evolution model stud- 
ied in 34.11 ). All models have Poisson error bars indicated. 
The four lower-z points at z = 4.25, 4.75, 5.25, and 6.75 con- 
tain N = 1, 3, 6, and 4 GRBs, respectively. The highest red- 
shift bin (z « 8.7 5) contains the two highest-redshift GRBs 
observed (z = 8.2. Tanvir et al.l l2009l: ISalvaterra et al.ll2009at 
z = 9.4, Cucchiara et al. 2011). 

Clearly the star formation rate densities estimated from the 
high redshift GRB rate through Equation [16] for physically- 
plausible models are considerably higher than those inferred 
from UV-selected galaxies. The model GRB rate to SFR ratio 
calculated in ^has a redshift-dependence that is constant at 
i}j ^ \ above z > 4, and the GRB-inferred /3*(z) at high red- 
shift is therefore large in this model. By consid ering the re- 
cently discovered GRBs at the highest redshifts dTanvir et alj 
I2009L ISalvaterra et alJl200"9ilCucchiara et al.ll201 ll) . we have 
extended this result to z ~ 9.5. The results are in contrast to re- 
cent clai ms for a low abundance of z > 10 star-forming galax- 
ies (e.g.. lBouwe ns et al. 2011; Oesch et al. 2011). To recon- 
cile /0*(z) estimates from both GRBs and UV-selected galaxies 
would require a dramatic evolution in ^(z). The physical ba- 
sis for such an evolution is unclear (c.f., ® . 

6. DISCUSSION: IMPLICATION FOR COSMIC REIONIZATION 

As the discovery of high-redshift galaxies (M cLure et al.l 
i2010; .Bouwens et al...2010b; .Oesch et al...2010.) and quasars 
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Figure 5. Implications of GRB-derived estimates for the high-redshift star formati on rate density, Panel a (upper left) shows p*{z) determined from UV-selected galaxies 

(gray points with error bars, see text) and coiresponding parametric SFR histories from'Robertson et aL'^2010^ (green hatched region). Also shown are values implied by the GRB rate 
assuming no evolution in '^{z) (red points), our model for GRB production in low-metallicity galaxies (blue points), and strong evolution in GRB production per unit star formation rate 
(^^(z) (X (1 +z)'^, orange points). The model points have been offset slightly in redshift for clarity and the model error bars reflect Poisson eiTors on the GRB rate in each redshift bin. 
If the GRB rate to SFR ratio evolves weakly beyond z > 4 (red and blue points), the rate of discovery of high redshift GRBs already implies a p*(z) much larger than that inferred from 
UV-selected galaxies. Evolution faster than "^(z) oc (1 +z)' would be needed to force agreement. Parameterized star formation histories consistent with the GRB-derived star formation 
histories in the constant 'I' and low metallicity star formation models are shown as black lines. With fiducial choices about the character of the stellar populations (Z ~ 0.2Zq), the 
dumpiness of the intergalactic medium (C = 3, upper line; C = 2.5 lower line), and the escape fraction of ionizing photons (/e^c = 0.06, upper line; f^^^ = 0.2, lower line) we can 
calculate the reionization history in panel b (upper right) implied by the GRB-derived high-redshift star formation rate (black lines) and compare with similar histories calculated 
by Robertson et al. i2010) determined from UV-selected galaxies (orange hatched region). The GRB-inferred star formation history would produce a large volume filling fraction of 
ionized gas extending to high-redshift. The path length through this ionized gas to the cosmic microwave background provides the optical depth to electron scattering in panel d 
(lower right). The ionization history computed from the GRB-derived star formation history would easily reach implied by the seven-year Wilkinson Microwave Anisotropy Prolye 
measurements i Komatsu et al. 201 1 ), and produce a much larger value than that similarly calculated from UV-selected galaxies (red hatched area; Robertson et al. 2010). While both 
the ionization history and the Thomson optical depth depend on specific model choices for /esc or C, the stellar mass density (panel c, lower left) is simply determined by the integral of 
the previous star formation rate density (panel a, upper left). The stellar mass density to z ~ 8 is shown as gray points with error bars I Gonz alez et al. 201 1), with the associated models 
by Robeitson et al. 1 2010, blue hatched region). The black lines in panel c show the stellar mass density implied by parameterizations of the GRB-derived star formation rate, which 
clearly exceed the stellar mass density at all redshifts. 

IO'^Mq yr"' Mpc~^. With metal poor stellar populations, a 
typical escape fraction of /esc ~ 0.3 and a clumping factor 
C ^ 2, the upper /0*(z) curve fully ionizes the intergalac- 
tic medium by z ~ 7 and recovers the Wilkinson Microwave 
Anisotropy Pr obe (WMAP) e lectron scattering optical mea- 
surement (e.g. jKomatsu et aTl bOll). 

Using the same form of Equation[T7]to parameterize the star 
formation rate density implied by the high-redshift GRB rate 
and a weak-to-moderate redshift-dependence of ^>{z) (Fig- 
ure |5] panel a, red and blue points), we find that parameter 



(iMortlock et aUlMTTl) reaches beyond z > 7, it is important 
to understand the potential role for star-forming galaxies in 
reioniz aton (for a review, see Robertson et al. 2010). For in- 
stance, [Robertson et al.l (gOlO) have calculated the reioniza- 
tion history, Thomson scattering optical depth, and stellar 
mass build-up for the star formation histories plotted in Fig- 
ure |5] (hatched regions). The star formation histories in these 
models have been parameterized using the formula 



a + b{z/cy 



(17) 



b = Q.IIMq yr" 



Mpc 

-3 



which 
model 



is the 'Robertson et al.' ('2010') generalization of the 
by Cole et al. (2001) to include a floor in the star 
formation rate. In Figure |5] panel a, the upper star for- 
mation history of the green hatched area has parameter val- 
ues a = 0.009Mq yr-i Mpc"\ b = O.IIMq yr"' Mpc"^ c = 
3.7, d = 7.4, and g = 10"^Mq yr"' Mpc"^ The lower star 
formation history model has c = 3.4, (/ = 8.3, and g = 



values of fl = O.OO7M0 yr"' Mpc 

c = 3.7, d = 6.4, / = 2.5 and g = 3 X 10"''Mo yr"' Mpc"- are 
representative as a "low-/ii(z)" GRB-derived model (lower 
black curve), while adopting d = 7.4 and ,g = [(4 x 10"^- 
10"^) X (z/3)+ 1O"-']M0 yr"' Mpc"-' provides representative 
"high-/3*(z)" GRB-derived model (upper black curve). 

With the same assumptions for the escape fraction, IGM 
clumping facto r , and stellar population model used by 
[Robertson et alJ (120101) . the star formation history implied by 
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the GRB rate for either a constant '^(z) or an evolution of 
the GRB rate to SFR ratio ^(z) that tracks star formation 
in low-metalUcity galaxies would fully reionize the universe 
by z ^ 14 and significantly over-predict the Thomson scat- 
tering optical depth (e.g., t 0.2). If instead high-redshift 
galaxies have metallicities Z ^ 0.2Zq with a Schaerer (20031) 
stellar population, the escape fraction is lower (/esc ~ 0.06), 
and the clumping factor is C ^ 3, the high- p^,(z) model im- 
plied by the GRB rate induces reionization by z ~ 7 and pro- 
duces a Thomson scattering optical depth which can match 
the WMAP value (upper black curve in Figure|5] panels b and 
d). For /esc = 0.2 and C = 3, the low- Pi,(z) GRB-derived model 
produces some what lower ionized gas volume filling factors 
and Thom son optical depths (l ower black curve in Figure |5]l. 
Models by Wvithe et al.l (120101) reach similar conclusions. 

Although the uncertainties associated with calculating jd*(z) 
from GRBs are large, an explanation of the possible dis- 
crepancy between the GRB-inferred /i*(z) and the abundance 
of high-redshift galaxies is warranted. The higher values 
deduced for the GRB-inferred /0*(z) compared to those in- 
ferred from UV-selected galaxies may strengthen the case for 
a steep luminosity function for the latter and hence closure 
on the hypothesis that intrinsically faint star-forming galaxies 
over 6< z < 12 were responsible for cosmic reionization (e.g., 
lYtiksel e t al."2008l: iK istler et al. l l20"09.) . However, as empha- 
sized by Robertson et al. (201^, there is a further constraint 
provided by the observed stellar mass density at z ~5-6. This 
constraint is determined from rest-frame optical fluxes de- 
duced from the Infrared Array Camera (IRAC) instrument on- 
board Spitzer Space Telescope. Importantly, the pAz) implied 
by the high redshift GRB rate appears unphysical in that it 
overproduces the observed stellar mass density af z > 5. Panel 
c of Figure |5] shows t he rest-frame optical s tellar mass den- 
sity determinations by'Gonzalez et al. (2011, gray points with 
error bars; s ee also Evies et al. 2005, Stark et al. 2007, 20091 
iLabbe et al. 2006, 2010land lGonzalez et al.l2010l) . These data 
may be compared with the integral of the parameterized star 
formation rate densities consistent by the high-redshift GRB- 
derived /d*(z) (Figure|5] panel c, black lines). While the escape 
fraction, clumping factor, and stellar population properties are 
not well known and may be sensibly varied within the broad 
uncertainties to match the observed Thomson optical depth to 
electron scattering, simultaneously accounting for the compa- 
rably low observed stellar mass density at high redshift is dif- 
ficult. While deeper IRAC observations and other techniques 
are improving these constraints (Gonzalez et al., in prepara- 
tion) it is important to recognize that these new estimates have 
confirmed that previo us efforts have p roperly accounted for 
incompleteness (e.g.. IStark et al.ll2007h and the stellar mass 
density measures are unlikely to increase substantially owing 
to future data. 

We thus conclude there is an important conflict between 
fairly reasonable assumptions about how the GRB rate 
maps to cosmic star formation and what we already under- 
stand about early star formation from UV-selected galax- 
ies. Can the two probes of early star formation be rec- 
onciled? Dust is an unlikely explanation given we would 
require all UV-selected star forming galaxies to be heavily 
extincted and most studies have utilized their UV contin- 
uum slopes to infer the absence of any significant reddening 
(lBou\ye ns et al. ll2010cl: iFinkelstein et aLllioiOt iMcLure et all 
l2011l: lDunlop et al.l2011h . The high-redshift UV galaxy lumi- 
nosity functions used to determine the star formation rate den- 
sity data in Figure |5] are taken from .Bouwens et al.. (.2010ah . 



At z = 8 the lBouwens et al.l (l2010al) luminosity function has a 
measured faint end slope of 7 = -2, and the luminosity den- 
sity at z = 8 shown in Figure |5] corresponds to a limiting ab- 
solute magnitude of Mab = -18. It has been suggested that 
extending the search to much fainter sources will help bridge 
the gap bet ween the galaxy - an d GRB-inferred P*{z) (se e Fig- 
ure 4 of K istler et al.ll2009l and lChoi & Nagaminell20ri ). but 
even with the very steep 7 = -2 the luminosity function would 
need to continue down to Mab = -(7 - 9) - slightly brighter 
than the globular cluster scale - for the UV-luminosity den- 
sity to recover the GRB-inferred star formation rate density. 
Deeper Hubble Space Telescope (HST) data will clarify the 
possibility. Adjusting the early stellar initial mass function 
(IMF) so as to increase the luminosity outp ut per un it pi^jz) 
would reduce the stellar mass density; for a ISalpeted ( Il955h 
IMF with an upper mass limit of M^ax = lOOM©, the min- 
imum mass of forming stars would have to increase from 
Mmia « O.IM0 to Mmin w 2.75 -4Mq to decrease the long- 
lived stellar mass by the required factor of f^, '^4-5. Con- 
ceivably, the IRAC fluxes at z ~5-6 are significantly contam- 
inated by nebular emission so that the deduced stellar masses 
are over-estimated. An adjustment to the IMF could provide 
as much as /^ times as much nebular emission contamination 
in the IRAC ban ds compared to th e nebular emission expected 
from a standard iSalpeteii ( 119551) stellar population. In this 
case, one factor of /* increase in the nebular emission contri- 
bution to the IRAC flux would arise from the increase in the 
number of Lyman continuum photons produced per unit star 
formation, and another factor of /, would arise from the de- 
crease in the contribution of long-lived stars to the rest-frame 
optical emission relative to a Salpete^ (119551) IMF. Luminos- 
ity evolution in the GRB population could also contribute by 
altering \l/(z) (e.g., [L lovd-Ronning et al. 2002; Firmani et al] 
2004; Kocevski & Liang 2006; Salvaterr a & Chincarinil2007F 
Salv aterra et al. 2009b; Virgili et al. 201 1 ), but this is also un - 
certain (e.g., (Butler et al. 2010; Wande rman & PiranI 12010). 
More likely, there are additional physical factors affecting 
the high redshift GRB production rate, i.e. beyond those 
of star formation and simple metallicity dependence consid- 
ered in this paper. Possible examples include the physics 
of angular momentum retention in GRB progenitors (e.g., 
iMacFadven et al.l 1200 Ih or the effects of the cosmic back- 
grou nd radiation tempe rature on the initial mass function 
(e.g.. lLarsonll998ll2005h . Regardless of what the cause might 
be, it is clear that the continued discovery and study of z >6 
GRBs promises to shed light on the reionization process and 
the assumptions made in the interpretation of star formation 
rates and stellar masses from UV-selected galaxies. 

7. SUMMARY 

Using the second Swift BAT catalog of GRBs 
(Sakamoto et al. 20 lib, the observations of d ark GRBs b y 
Perlev et all (120091) . iGreiner et af] (1201 Ih. an d Kruhler et al] 
go Hi) , and the GRB catalogs of iButler et a'l.. (J007, 2010), 
we have constructed the cumulative redshift distribution 
A^(< z\zmax) of 112 luminous (Liso > 10^' ergs s"') GRBs out 
to redshift z 9. Comparing with models of an evolving 
GRB rate to SFR ratio *(z), we find that A^(< z|zmax) 
constrains redshift evolution of the form ^'(z) oc (1+z)" to 
5, Q! 5, 1-5 out to Z max ~ 4. By extending the model of 
iKocevski et al.l (120091) to calculate ^'(z) from the evolution 
of the mass-metallicity relati on (|S avaglio et al. 2005), the 
star formation-mass relation (Drorv & Alvarez) i2008i) and 
the stellar mass function (.Drorv et al...2005.) . we find that the 
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presence of a host galaxy metallicity ceiling 12 + log[0/H]ciit 
above which GRBs are suppressed is highly consistent with 
the available data if 12 + log[0/H]crit < 8.85. The peak 
probability of consistency between the model and data occurs 
at l2 + log[0/ H1^iHt ^ 8.7, n ear the GRB metallicity ceiling 
suggested by . Modiaz e t al.l (12008). Using the method of 
lYuksel et all (120081) and Kistleretal. (2009j), we use the 



GRB rate at 1 < z < 4 (where the star formation rate density 
Pi,(z) is roughly constant) to estimate the at z > 4 

(including constraints from the highest redshift GRBs, e.g., 
Kawai et al. 2006; Tanvir et al. 2009; Salvaterra et al. 2009a; 
Cucchiara et al.ll201 1 ). We find that for constant to moderate 



(a < 1) redshift evolution in ^(z), the star formation rate 
density predicted by the observed high-redshift GRB rate is 
substantially higher than the observed fe(z) inferred from 
the abundance of UV-selecte d galaxies (e.g., iBouwens et alJ 
I2010bt iMcLure et al] |2010^, and would over-prod uce the 
obser v ed high-redshif t stella r mass density (e.g.. Stark et alj 
120091; iGonzalez et alj 1201 Ih . Rough agreement between 
the UV- and GRB -determined p*(z) can be achieved if the 
redshift-dependence of ^'(z) is as strong as a > 1 .5 at z > 4. 
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APPENDIX 

Our analysis relies heavily on the excellent GRB catalog s provided by lButler et al.l (12007b . iButler et alJ (120101). and the Secon d 
Swift Burst A l ert Tel escop e catalog ('Sakamot o et al.ll201 ik . Further, we make use of the important studies bv lPerlev" etaDdSOOl), 
iGreiner et alj ( 1201 Ih . and iKriihler et al. (201ll) that provided redshift constraints for a sample of dark GRBs. Below in Table 
[U we provide the union of these catalogs (162 GRBs in total) as a convenience to the reader only. The papers providing these 
catalogs should be consulted for important details on the fluence measurements and redshift determinations, along with further 
original works discovering and characterizing individual GRBs. The table presents redshifts, isotropic equivalent energies and 
luminosities, and burst durations as described in 32.11 Of these quantities, only the isotropic-equivalent energies for 29 GRBs 
(noted belo w) and the luminos i ties (f or all GRBs) are newly computed in this work. Most of the E^sn values reported (129) are 
taken from IButler et al.l (120071 l2010h . with 4 recent Ei^o values adopted from iSakamoto et al] (1201 ll noted below). Of the 29 
GRBs whose energies we comput e as Ej^o = A-ndjS, where di is the luminosity distance to a redshift z and S is the fluence, 21 
broad-band fluences a re taken from Butler et al.l(l2007ll20Toh and 8 10 - ISOkeV fluences (effectively lower limits) are taken from 
ISakamoto et al.l(l201 1). These GRBs are also noted in the table. Dark GRB redshifts (indicated as upper limits where appropriate 
and noted) are adopted from Perley et al. (2009), Greineretal. (2011), or Kriihler et al. (2011). Their reported energies and 
luminosities are calcu lated assuming the GRBs lie at the redshi ft upper lirn its. Burst durations mostly are taken from Table 1 of 
ISakamoto et al.l (1201 1,) . except for 3 values taken from .Butler et al.l (120071 12010. noted in table). For completeness reasons as 
our analysis makes use of luminous (Liso > 10^' ergs s~') GRBs (see ^2. Il l, and less luminous GRBs are therefore indicated with 
italicized names. 



Table 1 

GRB Catalog 



GRB 



£i,o [10=-^ ergs] rgo M Li,o [lO^'^ ergs .s-'] 



GRB 



£i» [10« ergs] 



'90 [s] 



[10" ergs s-'] 



050126 
050223 
050315 
050318 
050319 
050401 
050412 
050416A 
050505 
050525A 
050603 
050607 
0507 13A 
050730 
050801 
050802 
050803 
050814 
050820A 
050824 
050826 



1.29 
0.58 
1.95 
1.44 
3.24 
2.90 
<4.50'' 
0.65 
4.28 
0.61 
2.82 
<4.00'' 
<3.6()'' 
3.97 
1.38 
1.71 
0.42 
5.30 
2.61 
0.83 
0.30 



0.80-;- 
0.070t!;5« 
5.7+*-2 

4 6+^-^ 
32 0"^'^° 
101 



+18'3a 



o.io; 

16.01 

2 0+*'''' 
^■"-0.1 

50 0+^" ° 

J2 3+10.8a 

9 o;8 
1.8^j 
6>' 



20.0 
0.15 



0.03 



1+20.0 

m 

-on? 



48.0 
21.7 
95.6 
40.0 
152 
33.3 
26.5 
6.62 
58.9 
8.84 
22.0 
48.0 
94.9 
145 
19.4 
27.5 
88.1 
144 
240f 
24.8 
35.7 



"■''^'^-o.oio 
0051+" '"'" 

"■"'"^-0,0007 

n IR+<'19 

"■^°-0.003 

0.073+"^? 



0.13!j 
3.75!j 

21. i;- 

0.025 
1.43 



roll 

0.37+8:1 
8.68:6«5 

1 28 

7 54+§;ii 
' -2.32 
31+0 27 

o.osoit 

0039+" O"'? 
26+" " 
30+"^" 

0011+°-^'^^'^'' 



0.011+ 



070521 

070529 

070611 

070612A 

07072 IB 

070802 

07081 OA 

071003 

071010A 

071010B 

071011 

071020 

071025 

071031 

071117 

071122 

080129 

080210 

080310 

0803 19A 

0803 19B 



1.35''' 
2.50 
2.04 
0.62 
3.63 
2.45 
2.17 
1.60 
0.98 
0.95 

<5.00'' 
2.15 
5.20 
2.69 
1.33 
1.14 
4.35 
2.64 
2.43 

<2.2()'' 
0.94 



25.2" 



+22.0a 
8.8 

9.0!'" 



0.50;^; 
2.0; 
30.0; 
0.50; 
0.90: 

18.0 

2^4+18ia 

10.0+fpO 
428 



Mo 



if. 



3 Q+^-1 
1 Q+fl.8 



1+0.50 



0.30:^,0 

■^■^-0.9 
I+IO.S 

i^5a 
^-3.7 
400+2"" 



5.9; 
29.2 



38.6 
109 
13.2 
365 
337 
15.8 
9.04 
148 
6.32 
36.1 
80.9 
4.30 
241 
150 
6.07 
80.0 
50.2 
39.4 
352 
43.6 
125 



1 54+'"' 
29+0-29 

]2^'^^ 

"■^ -0.02 

0089+^' 008 
Hl+fi^f 

0^1 rHi! 



-+8:W5 



0.32; 
0.32 
0.041 

174^|f 

7 ^7+7.32 
-1,20 
11.0+]| 

0.096||i 

0.0080j|34 

0.85^;' 
47+0-« 

'^■^'-0.08 

057+^- '"2 
2i4*§f 

^■^ -0.27 
6 21+3 -10 
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Table 1 — Continued 



GRB 



£i„ [1052 ergs] Li,o [10« ergs s-'] GRB 



£i» [10« ergs] fgo [s] L,so [10« ergs 



050904 

050908 

050915A 

050922C 

051016B 

05 11 09 A 

0511098 

051111 

060108 

060110 

060115 

060116 

060124 

060202 

060206 

060210 

060218 

060223A 

060418 

060428B 

060502A 

0605 lOB 

060512 

060522 

060526 

060604 

060605 

060607A 
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"Eiso calculated from fluence provided by lButler et alJ (l2007ll2010h . 
''Dark GRB redshift limit from'Perl ev et al.1 (l2009h . 
'■Redshift from .Perlev et al. (.2009.) . 
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''Ei.o from Table 1 3 of lSakamoto et al.1 (1201 lb . 

''£'iso calculated from fluen ce provided by Table 2 of lSakamoto et al\ (1201 Ih . 
^Burst duration taken from 'Butler et al.1 12007ll2010l) . 
SRedshift from Greiner et al. (201 1|. 
''Redshift from Kruhler et al. (201 li 



'Dark GRB redshift limit from. 




